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Abstract 

We present several different ways of describing a spherically symmetric spacetime and 
the resulting metrics. We then focus our discussion on an especially useful form of the 
metric of a spherically symmetric spacetime in polar-areal coordinates and its properties. In 
particular, we show how the metric component functions chosen are extremely compatible 
with notions in Newtonian mechanics. We also show the monotonicity of the Hawking 
mass in these coordinates. Finally, we discuss how these coordinates and the metric can 
be used to solve the spherically symmetric Einstein-Klein-Gordon equations and how this 
will be useful in our future work. 



1 Introduction 

Spherically symmetric spacetimes are an important case in the study of general relativity for 
a number of reasons. Foremost among them is that it is often a good starting point in the 
study of a problem in general relativity. For example, one of the first projects undergone 
in general relativity was to compute nontrivial spherically symmetric spacetimes that are exact 
solutions of the Einstein equation. This resulted in the discovery of the Schwarzschild spacetime, 
which is by far the most important spherically symmetric solution to date, and later Birkhoff 's 



theorem about Ricci flat or vacuum spherically symmetric spacetimes. 14 Spherically symmetric 
spacetimes also create a situation where the dynamics of the system are less complicated by 
effectively reducing a 4-dimensional solution to a 2-dimensional one. This accessibility makes 
using spherically symmetric spacetimes all the more attractive as a starting point. Finally, while 
Birkhoff's theorem classifies all vacuum spherically symmetric spacetimes, there are still some 
nonvacuum spherically symmetric spacetimes that are interesting as well, both from a physical 
and mathematical standpoint, such as Dwarf spheroidal galaxies, which are dominated by their 



Dark Matter halos and also closely approximated by spherical symmetry 13 , and spherically 
symmetric scalar fields. 

Since spherically symmetric spacetimes are still of interest, it would be useful to find an 
efficient way to describe and model these spacetimes in a general sense. In particular, it would 
be useful to have a form of a general spherically symmetric metric that is well-suited to numerical 
evolutions of the Einstein equation. There are many possible forms of the metric of a general 
spherically symmetric spacetime to use and the purpose of this article is to first collect these 
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metrics, discuss tlie advantages and disadvantages of using each one, and then present a metric 
that is not only extremely well-suited to numerical evolutions, but also is described in terms 
that are very natural analogues to the low field or Newtonian limit. Most of these metrics can 
be described well within the established framework of numerical relativity and as such, we will 
use that framework in the following. 



2 Metrics for a Spherically Symmetric Spacetime 



The study of numerical relativity is devoted to devising ways of evolving the Einstein equation 
in order to solve for the components of the spacetime metric in different situations of interest 
as well as actually conducting such numerical experiments and comparing them to real data. 
This evolution usually takes place in a spacetime which is foliated hj t = constant spacelike 
hypersurfaces. The common method in numerical relativity is to decouple the time component 



from the space components into what is commonly called the (3+l)-formalism. [Tip, 10 



The framework for this formalism is, as stated before, a spacetime, A^, foliated hy t = 
constant spacelike hypersurfaces described by a Riemannian 3-metric 7 which may change with 
time. Note that such a foliation is possible for any globally hyperbolic spacetime. (l, 14 Consider 
a coordinate chart on U N, {t, x^, x^, x^} where the dt is timelike and the d^j are all spacelike. 
Now consider an observer starting on the t = to hypersurface at the coordinate (to, 2^0)- This 
observer then travels to another infinitesimally close hypersurface t = to + dt to the coordinate 
(to + dt,xl + dx) as in Figure nj . 



to + dt 



to 



(to + dt, 4 - f3'dt) 



(to + dt, X* 




Figure 1: Inifinitesimal distance in a t = constant hypersurface foliated spacetime. 

The observer has now traveled an infinitesimal distance ds. We can measure the square of this 
infinitesimal distance, ds"^, using the analogue of Pythagorean's theorem and this will give us 
the line element form of the metric. If another observer travels normal to the hypersurface from 
(to,Xo) to the hypersurface t = tg + dt, since the normal direction is not necessarily the same 
direction as the t coordinate direction, it will arrive at the coordinate {tQ + dt,XQ — (3^dt). We call 
the 3-vector field, /3, the shift vector because it measures the spacelike shift of the coordinates 
while traveling normally. Note that the components of /3 can vary with all the coordinates. This 
normal observer, having traveled in a timelike direction, has experienced some proper time dr 
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which is some multiple of the change in time coordinate, that is, 

dr = adt. (1) 

Hence the length of its normal movement from one surface to the other is adt. The value of a 
can vary with all of the coordinates making it a function on the manifold, which is called the 
lapse function since it measures the lapse in proper time compared to coordinate time. To get 
the length in the spatial direction between where the normal observer ended up and where the 
original observer did, we need only to use the metric on the hypersurfaces and the difference of 
the two space coordinates. This difference, for each z, takes the form 

{x{ + dx^) - {xl - P^dt) = dx^ + P^dt (2) 

and so the length squared of the spatial movement will be 

'yjkidx^ + P^dt){dx^ + P^dt) (3) 

where we have implemented the Einstein summation convention. Then using the generalized 
version of Pythagorean's theorem and recalling that t is a timelike direction we get that 

ds^ = -a^ df + -fjkidx^ + P^dt){dx'' + P'^dt), (4) 

which is the line element of the metric. This is the most general form of the metric for any 
foliated spacetime, that is, all metrics of a foliated spacetime can be written in this form. Note 
that we will often write g instead of ds^ to refer interchangeably to the metric and the line 
element. 

From this metric, and a choice of slicing condition, a complete system of partial differential 
equations that evolve the Einstein equation can be constructed. This system is often referred to 
as the ADM formulation of general relativity in reference to the authors of the paper in which 
it was first introduced. |4] It involves evolution equations of both the metric and the extrinsic 



curvature of the t = constant hypersurfaces. |4||6, 10 These equations are very commonly used 
in numerical relativity, but we find they overcomplicate the situation in spherical symmetry, 
which is why we have elected not to use this formulation of general relativity directly. 

If we know more about the spacetime in question, we will be able to determine more of 
the components of the metric. In a spherically symmetric spacetime, with coordinates, r, 6, if, 
chosen so that 6 and if are the polar-angular coordinates on the hypersurface, the shift vector 
must be completely radial so that the metric remains invariant under rotations. In this case, we 
will denote the radial component of the shift vector by simply /3. Moreover, the 3-metric can be 
written as 

J = 'jrr dr"^ + 'yee da^ . (5) 

where da^ = dO"^ + sin^ 9 dip^ is the standard metric on the unit sphere. This implies that the 
most general spherically symmetric metric can be written in the form 

g = —o? dt^ + ^rridr + /3 dt^ [dr + /3 lit) + do^ 
= — (a^ — 'jrrP'^) dt^ + -yrrPidr dt + dt dr) + •jrr dr"^ + 76)e da'^. (6) 
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Note that all the metric component functions can only depend on t and r due to spherical 
symmetry. For convenience, we will define two positive functions a{t,r) and q{t,r) so that 

a(t,r)^ = 7rr and q{t,r)^ = jeg. (7) 

Then we can rewrite ^ as 

where a,a,P,q are functions of only t and r. It is important to note here that since we have 
not defined the coordinates t and r geometrically yet, there remain two degrees of freedom left 
in this metric. There are several different choices that can be made in this regard. We mention 
the most common here, but there is a very useful and more extensive list of several choices that 
can be made in a more general setting in Gourgoulhon's recent book, [l] 

There are three rather common slicing conditions that are often used in many settings, all 
of which place a condition on the lapse function a. These conditions are the maximal slicing, 
harmonic slicing, and geodesic slicing conditions. 

Under the maximal slicing condition, one requires that each t = constant hypersurface be a 
maximal hypersurface, that is, it has zero mean curvature, 

H = 0. (9) 

The metric stays of the form in equation ([6]), but since the evolution equations in the ADM 
formulation evolve the components of the second fundamental form and H is the trace of the 
second fundamental form, this places a constraint on some of the evolution variables, which can 



be used to simplify the evolution equations. [l|[5, 11 Note here that after making this choice, 
there remains one more degree of freedom which can be used to constrain the r coordinate. 

Harmonic slicing requires that the coordinate function t be a harmonic function under the 
metric g. That is, Ogt = 0, where Og is the dAlembertian or Laplacian operator with respect 
to the metric g. This is often accompanied with the condition that the hypersurfaces remain 
orthogonal to the time direction, which uses the remaining degree of freedom, and indeed some 
refer to both of these choices together as harmonic slicing. In the case that both conditions are 
satisfied, this would yield a metric of the form, 

g = -a^ dt^ + dr"^ + da"^ (10) 

with the added condition Dgt = 0, which can be used to compute an evolution equation for the 
lapse function a.j^ 

Geodesic slicing requires that movement along the curve ^ = (t, 0, 0, 0), which is given in our 
coordinates, be geodesic, that is, coordinate observer worldlines are geodesies. This requirement 
is satisfied by choosing 

a = constant and 13 = 0. (11) 

However, the most reasonable choice for the constant is 1, since the condition a = 1 on the lapse 
function has the added implication that normal observers proper time is the same as coordinate 
time and in fact, since (3 = 0, normal observers are coordinate observers. (T|[3||5] This results in 
a metric of the form 

g = -dt^ + dr'^ + q^ da^. (12) 
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While this choice seems very attractive at first, since it either eliminates or greatly simplifies the 
evolution equations, it does have a tendency to develop coordinate singularities when evolved 
in time and hence must be used with caution. [5] Note also that this requirement is sometimes 
reduced to simply a = 1 without necessarily requiring the shift parameter or vector to vanish. 

We can alternatively use the degrees of freedom to make choices concerning the r-coordinate. 
We will mention three such choices here that are standard in the study of spherically symmetric 
spacetimes. 

The first is the choice of normal slicing. That is, choose the t coordinate so that the dt vector 
field is always normal to the hyper surf aces. This choice makes all normal observers coordinate 
observers as well. This is equivalent to choosing the shift parameter or vector to be identically 
and results in a metric of the form 

g = -^2 ^^2 ^ ^2 ^^2 ^ ^2 ^^2_ ^^3^) 

This choice has already been mentioned above as it is often coupled with harmonic or geodesic 
slicing conditions. |T| Since the coordinate vector fields on the hypersurfaces were already or- 
thogonal, this results, as seen above, in a diagonal metric. 

The next choice is to require that the metric on the hypersurfaces to be conformal to the 
flat metric. This amounts to choosing the function g in (|8]) to satisfy q = ra, which would make 
the metric become 

5, = - (a^ - a^/?^) dt"^ + a^P{dr dt + dt dr) + [dr"^ + rfa^) . (14) 

This choice is referred to as isotropic coordinates. It is often coupled with the normal slicing 
choice above, which uses both of the degrees of freedom and results in a metric of the form 

g = -a^ dt"^ + {dr^ + da^) . (15) 

A common example of the normal-isotropic case is the Schwarzschild metric in isotropic coordi- 
nates. (l4] 

The last choice we mention in the (3+l)-framework is to give the coordinate r geometric 
significance by choosing it to be the areal coordinate. That is, choose the coordinate r so that 
the area of each metric 2-sphere on the hypersurface is exactly Airr"^. This choice requires that 
q = r. Additionally, it is almost always accompanied with the normal slicing choice above, 
again using both degrees of freedom. This results in the polar-areal coordinates on a spherically 
symmetric spacetime and yields a metric of the form 

g = -a'^ dt'^ + a'^ dr^ + da^. (16) 

This chart is probably the most familiar form of a general spherically symmetric metric, and is 
the chart most commonly used when introducing the Schwarzschild spacetime. And rightly so, 
as it has some very clear advantages. For one, the r-coordinate 's role is analogous to its role 
in a fiat spacetime. Additionally, these coordinates give the Einstein curvature tensor a very 
simple form. However, they may not be well suited to dealing with high gravitational fields as 
we would likely run into the same limitations that the polar areal metric has in describing the 
Schwarzschild spacetime inside the Schwarzschild radius. This is not much of a problem for us, 
since we are mostly interested in the low field limit. 
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Finally, we present one more useful coordinate system and metric for a general spherically 
symmetric spacetime that does not fit into the (3+l)-formalism framework, nor does it depend 
on the ability to foliate the spacetime, but is useful from a theoretical standpoint nonetheless. 
In these coordinates, we still choose to use the polar-angular coordinates 9 and 99 to describe 
the rotations, but instead of separating the time and radial coordinates, we choose coordinates 
u and V such that and are nonparallel future pointing null vectors. This coordinate system 
is descriptively called null coordinates. In this case, the metric on the spacetime takes the form, 

g = A (du dv + dv du) + da"^ (17) 

where A and Q are functions of only u and v. Note that there are no extra degrees of freedom 
with this metric as all coordinates have been well-defined. While these coordinates are not very 
well suited to numerical evolutions, they can be very useful for theoretical discussions about the 
spherically symmetric spacetime. For example, in these coordinates, it is very straightforward to 
prove the monotonicity of the Hawking mass given the dominant energy condition. Additionally, 
it seems to perform well when in the presence of high gravitational fields. In the Schwarzschild 
case, these coordinates are known as the Kruskal coordinate system and is the system generally 



used to describe the region inside the Schwarzschild radius. 14 



All of these different coordinate choices and different forms of the metric on a spherically 
symmetric spacetime have advantages to them. However, for the problem of numerically evolving 
a spacetime metric in a low gravitational field, we find that the polar-areal coordinates are 
the best suited due to the very simple system one gets from the Einstein equation. As such, 
polar-areal coordinates is the coordinate system that we will use throughout this paper, but, 
in addition, we will introduce new variables that will give the metric a different form. This 
new form of the metric will result in the added advantages that the new metric functions have 
very clear analogues in the Newtonian or low-field limit and the Einstein curvature tensor will 
become even more simplified. 



3 A Newtonian-Compatible Metric 

Let A be a spherically symmetric, (3+l)-dimensional, asymptotically Schwarzschild (that is, it is 
Schwarzschild as r — )■ 00) spacetime with a Lorentzian metric g and be foliated by t = constant 
spacelike hypersurfaces. Choose the polar-areal coordinate system discussed above so that the 
t-coordinate direction is always normal to the hypersurfaces, the radial coordinate r is the areal 
coordinate (that is, so that the metric sphere of metric radius r always has a surface area 
of 47rr^), and 9 and are the usual polar-angular coordinates on the hypersurface. In these 
coordinates, the metric g has the line element form 

g = -a{t, rf dt^ + a(t, rf dr'^ + da'^ (18) 



Now we define the functions V{t,r) and M(t,r) as follows, 

r f a? — I 
2 



V = \oga M=-(— ^) (19) 



and rewrite the metric (18) in terms of these functions to obtain. 



g = -^^v de+[l-—\ ' dr' + r'da' (20) 
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Note that we will always be interested in the low-field limit, where M <^ r, which necessarily 
requires that M{t, 0) = for all t. 



3.1 Compatibility with Newtonian Physics 

Here, we wish to compute several important properties about this metric and the physical 
interpretation of both V and M. To physically interpret this metric, we will introduce the 
Einstein equation, but first, we present a property that doesn't require the Einstein equation, 
the proof of which can be found in Appendix |X} 

Proposition 3.1. The function M{t,r) is the spacetime Hawking mass of the metric sphere, 
St,r; for any given t and r. 

This suggests that we should interpret M as the mass of the system. However, there is 
an even stronger reason to do so, which we will investigate later. In order to prepare for that 
discussion, we need some preliminary results first about the relationship between the metric and 
its stress-energy tensor. To begin, consider this metric as a solution to the Einstein equation 

G = SttT. (21) 

To facilitate our discussion, we will compute the Einstein curvature tensor of this metric in 
certain directions. To that end, define the following unit vector fields. 



ut = e^dt Ur = \ 1 - — Or i^e = - de Vu, = — ^ <9^ (22) 

V r r rsmu 

Note that at every point, p E N, except the coordinate singularities r = and 9 = ±n (that is, 
all points where all the vector fields are well defined), these vector fields form an orthonormal 
basis of TpN and hence are a frame field. The Einstein curvature tensor is defined as 

G = Ric-^Rg (23) 

where Ric and R are the Ricci curvature tensor and the scalar curvature of the spacetime 
respectively. Since both the Ricci curvature tensor and its trace R are present in this equation, 
we will need to know a few results about the Ricci curvature in these coordinates. We have 
the following lemma and subsequent corollary. While the proof of the corollary is short and is 
presented here, the proof of the lemma can be found in Appendix [A} 

Lemma 3.2. The only nonzero components of the Ricci curvature tensor in the basis and 
the scalar curvature are as follows. 

1. Ric(z/i,z/i) = ( v„ + v;2 + — f 1-— j + ^ f y-M, 

VtMt-Mtt L 2My^ 3M2 2M 
2M, / 2M^ "^/^ 



2,. Ric(z/(, Vr 
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3. Ric(z/„ Ur) = -{Vrr + K ) (^1 " 

VtMt - Mtt 



2M 



re 



2V 



2M 



ly 2 ^2X^ 



Vr 



1 

2M 

r 



Mr 



2M 



-2 



5. = -2 Vrr + V; + 



r2 . 2K 



1 - 



2M\ 2Vr fM 



{—-Mr\ + 



2{VtMt - Mt, 



re 



2V 



1 - 



2M 



-1 



+ 



6M2 



1 - 



2M 



Corollary 3.3. The only nonzero components of the Einstein curvature tensor in the Vr) basis 
are as follows. 



1- G{yuvt)-. 

2. G{yuVr) 

3. G{Vr,Vr) 



2Mr 



2M 



-1/2 



2M 2Vr 
■— ^ + — 



1 - 



2M 



2M 



fM 



+ - 



Mr 



\ r 



+ 



VtMt - Mt, 



re 



2V 



1 - 



2M 



ml 

^ 2g2 V 



1 Vr 



2M 



Proof. Since {vf, i^r, T^e, ^tp] form an orthonormal basis everywliere. We liave tlie following. 



G{ut 


^t) 


= Ric{vt, 


^t) 


- ^diyu vt) -- 


= Ric{vt,vt) + 


R 

~2 


G{vt, 


Vr) 


= Ric{vt, 


Vr) 


- -^di^t, Vr) 


= Ric{vt, Vr) 




G{yr, 


Vr) 


= Ric{vr, 


Vr) 


- ^9{k,k) 


= Ric(z/r, Vr) - 


R 

~ ~2 


G{ve, 


I'e) 


= Ric{vg, 


I'e) 


- ^gi^e, ve) 


= Ric{vg, Vg) - 


R 
" ~2 


G{iyp, 


v^) 


= Ric{v^ 




R , 


) = Ric{v^, v^) 


R 

~ ~2 






= Ric{vg, 


ye) 


- f = G{ve. 


ve) 





(24) 
(25) 
(26) 
(27) 

(28) 



We then use Lemma |3.2| to substitute in the values for Ric and R. The rest is algebra. These 

sis 

□ 



are the only nonzero components by Lemma 3.2, equation (23), and the fact that the z/^ basis 
is orthonormal. 
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Next, we define the function fi{t,r) to be the energy density of an observer ai p = it,r) 
moving through the shces with 4- velocity Uf. That is, 

fi = T{ut,ut). (29) 

We now have enough information to prove the following proposition, which contains the promised 
stronger reason for interpreting M as the mass inside each metric sphere. 

Proposition 3.4. For a fixed t and r, M{t, r) is the flat integral of the energy density, fi{t, r), 
over the ball Ef r of radius r at time t. 



Proof. By the Einstein equation (21), equation (^29|, and Corollary 3.3, we have that 



G{vt,vt) = 8nT{iyt,iyt) 
,2 ^'^f^ 



2Mr 



r 

„2 



Mr = 47rr> (30) 

Since L dA = 47rr^, where Sj,, is the sphere of radius r at time t, we have then that for a 
fixed t and r, 

M{t,r)= / ATTs'^fi{t,s)ds = / / fi{t,s)dAds (31) 

Jo Jo JT,t,s 

l*T /»27r /"TT 

= 11 / fiit,s)s^ sine de dip ds (32) 
Jo Jo Jo 

/i(t, s) dVo, (33) 



where dVo = s"^ sine de dip ds is the fiat volume form on the ball of radius s, and we have 
introduced s as a dummy integrating variable in the "r" position. Thus M{t,r) is the fiat 
volume integral of the energy density over the ball of radius r. □ 

Note that M is not the integral of rhw energy density with respect to the metric's volume 
form, dV = (1 — 2M / s)~^^'^ s'^ sine ds de dip, but rather the following is true. 



M(t,r)= / i^{t,s)\ll-^-^^dV. (34) 

However, in the Newtonian limit, M ^ r, the above integral becomes approximately the integral 
of the energy density with respect to the metric's volume form over the ball Et^r of radius r. 
Thus referring to M{t, r) as the mass inside the metric sphere of radius r at time t not only 
makes sense from a geometrical point of view given the Hawking mass, but also from a physical 
point of view. 

Next we define the function P to be the pressure in the stress-energy tensor for an observer 
at (t, r), that is, let 

P = T{l^r,l^r)- (35) 



Then we can use Corollary 3.3 and the Einstein equation (21 ) to prove the following proposition. 
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Proposition 3.5. In the Newtonian limit, where P = and M <^ r, we have that AV = 47r/i, 
where A is the flat Laplacian on M^. 



Proof. By the Einstein equation (21), equation (^35|, and Corollary 3.3 we have that 



2M 2Vr 

■— ^ + — 





S7iT( 




2M\ 






= 


8nP 




r J 








0- 




r 








2M^ 


r'^Vr = 




r 



-1 



2M 



+ 87rP 



(M + Anr^P) . 



(36) 



Note that in the Newtonian limit, where P = and M ^ r, this equation is approximated by 

M. (37) 



Also, the metric on the hypersurface, under these assumptions, is approximately the polar-areal 
metric on M^. Note that the Laplacian on of a spherically symmetric function / is given by 



A/ 



1 d 



d^f 2df _ 1 d 



dr 



(38) 



Then applying the operator to (37) and using equation (30) yields 



J.2 Qy 



AV 



dr 

4m. 



(M) 



AV = Att/j, 



(39) 

□ 



Equation (39) is Poisson's equation and is the defining equation of the gravitational potential 
in Newtonian mechanics. Moreover, equation (37) reduces to Vr = M/r"^ which yields the inverse 
square law for Newtonian gravity. Then we can interpret V as the analogue in our scenario of 
the Newtonian potential. The interpretation of M and V via Propositions 3.1, 3.4, and 3.5 are 
what we mean by saying that the metric (20) is Newtonian compatible. 



3.2 Other Useful Properties of this Metric 

In this section, we produce two additional useful results which are readily obtainable with these 
coordinates and metric functions. The first is the well-known result of the monotonicity of the 
Hawking mass in spherical symmetry, which is made particularly straightforward using this form 
of the metric. It follows almost immediately from Corollary 3.3[ 
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Proposition 3.6. // the spacetime satisfies the dominant energy condition that G{X, Y) > for 
every future-pointing causal vector fields, X and Y , then, whenever 2M{t,r) < r, the Hawking 
Mass, M{t,r), is monotonically increasing in any non-timelike direction for which the radial 
coordinate increases. 

Proof. The vector field ut is the future-pointing timelike unit vector field in the t direction and 
the vector fields i>r + ut (future-pointing) and z/^ — ut (past-pointing) are null vector fields in the 
null directions where the r coordinate increases. By Corollary |3.3[ we have that 

/ 2M\^''^ Mf 
{u^ + u,){M) = Mr[l~—j 

_ r2 / 2My/^ / 2Mr 2Mt f 2M^ 



rW 2M^ 



2 V r 



rW 2M^ 



2 V r 



GiUt, Vt + Vr) (40) 



The last factor here is positive by the dominant energy condition since both Vt and Vi + are 
future-pointing causal vector fields. Since > and 2M < r everywhere, it must be that 
(z/j. + z/f)(M) > everywhere as well. By the same corollary, we also have 

(i/,-z/t)(M) = Mr (^1 

_ r2 / 2MV'^ (2Mr 2Mt ( 2M^'^'^ 

f 2M^ 




2 \}-—] {G{ut,ut)-G{ut,ur)) 



f 2M^ 



2 V r 



G{vt,vt-Vr) (41) 



The last factor here is positive by the dominant energy condition since vt — is also future- 
pointing causal. Then, as before, it must be that (z/^ — z/t)(M) > everywhere. Since both 
(z/j. + z/t)(M) and (z/^ — z/t)(M) are both nonnegative everywhere, any positive linear combination 
of the two is also nonnegative, which is the desired result. □ 

This next property will be useful later when we consider the Einstein- Klein- Gordon system. 
By Corollary |3.3 and the Einstein equation, any stress-energy tensor T for this spacetime could 



have no other nonzero components than the corresponding components of the Einstein curvature 



tensor presented in Corollary 3.3 Thus we define the following functions, two of which have 



already been introduced in equations (29) and (35), 



/x(t, r) = r(z/j, z/t) p(t, r) = T(z/t, v^) (42a) 

P(t, r) = r(z/„ z/,) Q(t, r) = T(z/,, z/,) = T(z/^, v^) (42b) 
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and these functions account for all the nonzero components of T in terms of the orthonormal 
frame {z/^, i/,., z/g, i/<^}. Then we have the following result which follows directly from the fact 
that Vg ■ T = 0, as with all stress-energy tensors, but whose lengthy proof we will leave to the 
appendix. 

Proposition 3.7. // the components of the stress-energy tensor T are known for all t and r 

and the functions M{t, 0) and V{t, 0) are known for all t, then solving the resulting initial value 



problem on the ODE's obtained from (30) and (36), namely 



Mr = 47rrV 



Vr 



1 - 



2M 



M 



+ AnrP 



) 



(43) 
(44) 



for each value oft, will solve the entire Einstein equation. 



4 The Einstein-Klein-Gordon Equations 

A good example of the usefulness of this choice of metric is in the evolution of the Einstein- 
Klein-Gordon equations, which is the Einstein equation in the presence of a nontrivial scalar 
field. In this paper, we will work specifically with a complex valued scalar field. 

Consider a spherically symmetric spacetime N with Lorentzian metric g as described in the 
previous section. Let / : — )■ C be a spherically symmetric complex valued scalar field. Then 
the Einstein-Klein Gordon equations in this case are 

( df®df + df®df (\df\\^.^2\ \ .... 
G = Snfio I — I — + 1/1 j gj (45) 

= TV (46) 

where /iq is simply some constant that controls the scale of the system and is not to be confused 
with the energy density fi{t, 0) at the central value. The value of /iq is unimportant to the 
qualitative behavior of the solutions and can be absorbed entirely into /, if desired. By contrast, 
the parameter T is a fundamental constant to the equations upon which the qualitative behavior 
of the solutions depend. The real valued version of these equations, that is, where / is a real 



valued scalar field, is derived in a paper by Bray. [9] Equations (45) and (46) are the Euler- 
Lagrange equations of the action involving the scalar field given in Bray's paper if the scalar field 
is assumed to be complex valued instead of real valued (this follows the Lagrangian formulation 
of general relativity given in an appendix in Wald's book[l^). One can readily see that we 
recover the equations in Bray's paper, if we assume that / is real valued, which implies that 



f = f. Note also that equations (21) and (45) imply that for a complex scalar field 



, df®df + df®df (\df\\...2. , 
T = /io I I ^ + 1/1 \9]- (47) 

With these equations, we will construct a system of partial differential equations which can be 
used to numerically evolve the scalar field / and the metric from consistent initial data. Since 
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in the previous section we have already computed the components of the Einstein curvature 
tensor, in order to construct such a system of PDEs, we need to compute the components of 
the stress energy tensor and the Laplacian in the metric g. Before we do, we will define the 
function, r), by the equation 

2M^ "^/^ 



1 - 



(45 



This is done to make the resulting system of PDEs first order in time and results in a more 
convenient choice than choosing only p = ft- We now have the following lemmas. 

Lemma 4.1. The following are true. 

2M\ \frf + H 



1. r(z/„z/i) = /"o| l/r+(i 

2. T{ut,Ur) 



^2 



3. T{iyr,i'r) 



4. T[ve, ve) = T[v^, v^) = -/iQ |/| + 1 



The proof of this lemma is left to the appendix. 



Lemma 4.2. The components of the Einstein equation in the u^j directions result in the following 
PDEs. 



87rrVoe^ / _ 2M\ 
I r J 



2M\ \frf+\pf' 



r 

3/2 



^2 



Re{frP) 



VtMt 



1 



2M 



-re 



2V 



M 



1 - 



2M\ iM'+br 



2M 



r 

iM 

\ r 



(49) 

(50) 
(51) 



1 Vr 
+ — 

r 



2M 



3M2 



2M 
r 



iTTjJo 1 - 



2M 



\ff+ I 



2M\ \ fr\'~\pf 
'^2 



(52) 



Proof. This follows directly from the Einstein equation (21), Corollary 3.3, and Lemma 4.1 and 



the algebra necessary to solve for each of the quantities above. Equation (49) follows from the 
z/(, Vt component of the Einstein equation, (50) follows from the Vf, component, (51) from the 
z/j., Pr component, and (52) from the pq, uq or i/^, component since they're identical. □ 



13 



Alan R. Parry 



A Survey of Spherically Symmetric Spacetimes 



The Klein-Gordon equation and the fact that 



9\ 

yield the following lemma. 

Lemma 4.3. The Klein- Gordon equation in this metric yields the following PDF. 

2M^ "^/^ 



r V r \ V r 



Pt = e' -r^f 1 



Proof. To compute the Laplacian, we'll first need the quantity y\g\- We have 



1^1 



Then since / = f{t,r) and g is diagonal, we have by equation (53) 

-1/2 



9 gjj^ Q 



r sin 6 g ft 



+ a e^ 1 



2M 



-1/2 



sm9 g^^ fr 



1 



2M 



k -e 



2M 



-1/2 



sin ^ 



+ dr\e''\ll-—r'sm9fr 



2M . . 2 

1 I Ot [—pr sin 6'j 



2A/f / 2M 
+ 2re^ W 1 sin ^ + sin ^ 9, e^/, W 1 



-Pie \ 1 



2M 

+ ^11 



2M 2fr 



+ e 



2M 



a. |e^/,Wl 



2M 



Then the Klein-Gordon equation (46) becomes 



-Pte 



/I _ !M + fi _ m) -^JLa^ e'7.,/i I - 

f T 1 V 



r 



Solving for yields the desired result. 



(53) 



(54) 



(55) 



(56) 



TV (57) 



□ 
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Then rewrite (48) as 



ft =pe \ I 



2M 



{5t 



to yield an evolution equation for /. 

We now include a lemma that shows that the T given in (47) has the property of all stress- 
energy tensors that Vg ■ T = 0. 



Lemma 4.4. Equation (46) implies that Vg ■ T = 0, where T is given by (4'^). 

This lemma's proof is also included in the appendix. Now we can prove the following. 



Lemma 4.5. If f,p^ M^V solve (49), (51), (54), and (58), then they also solve (50) and (52). 



Proof. This follows from equation (42), Lemmas 4.1, 4.2, 4.4, Proposition 3.7, and the fact that 
we can use (54) and (58) to evolve the components of the stress-energy tensor through time as 
long as we use (49) and (51) to compute the metric components at each new time. If we do 
these two procedures simultaneously, it is effectively the same as knowing the components of 
the stress-energy tensor for all values of t, which is required by the hypothesis of Proposition 
[321 □ 



Lemmas 4.2 4.3 and 4.5 and equation (58) prove the following proposition. 

Proposition 4.6. Solving the following PDEs for f{t,r), p{t,r), V{t,r), and M{t,r) with 
consistent initial data is necessary and sufficient to solve the Einstein-Klein- Gordon equations 
(45) and (46) with a metric of the form on the same initial data. 



Mr 



Vr 



ft 



Pt 




2M\ \fr? + \p[ 



^2 



-47rr/xo I/I'- 1- 



2M\ |M' + |p 



^2 




-1/2 



+ 



Vr 



1 - 



2M 



+ dr eVr\/l- 



2M 



(59a) 
(59b) 

(59c) 
(59d) 



5 Future Work 

Of all the different coordinate systems and metric functions that we have used, this Newtonian 
compatible choice yields by far the simplest system of PDEs that we have come across that 
solve the spherically symmetric Einstein-Klein-Gordon equations. Solutions to the Einstein- 
Klein-Gordon equations could possibly yield a geometric description of dark matter. This form 
of dark matter, which we call wave dark matter, has been studied for years, but has recently 
been well motivated geometrically by Bray. ^ In much of our future work, we plan to study the 
geometric and gravitational properties of solutions to the Einstein-Klein- Gordon equations in 
order to compare them to observations. 
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One of the best places to start is with dwarf spheroidal galaxies, which are small, roughly 
spherically symmetric, and dark matter dominated. The system (59) which solves the spherically 
symmetric Einstein-Klein- Gordon equations will be particularly useful in this regard. We will use 
this system to study both static and non-static solutions to the Einstein-Klein- Gordon equations 
and compare their predictions to observations of the velocities of the stars in these galaxies. To 
this end, we have created Matlab code to take initial conditions of the scalar field /(O, r) = /o('") 
and p{0,r) = po{r) and evolve the entire Einstein- Klein- Gordon system indefinitely. The code 
works on relatively simple finite difference techniques and will be discussed in another work. 
The results of these endeavors will help begin to answer the question of whether or not wave 
dark matter correctly describes the gravitational effects of dark matter and will also help us 
determine the value of T, the fundamental constant in the Einstein-Klein- Gordon equations. 
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A Proofs 

In this appendix, we present the proofs which were omitted in the above sections. Every proof of 
a proposition or lemma from above will refer back to the number of the statement it is proving. 
We also introduce a new lemma here which is useful for proving some of the above statements. 
We will present all of the omitted proofs here in the order with which their corresponding 
statements are given above. 



Proof of Proposition 3. 1 The Hawking mass of a metric sphere is defined as 



m^,(Si,,,) = A/^ (l-T^ / g{ll,E]dA\ (60) 



167r \ IGtt 



where \T,t^r\ is the surface area of the metric sphere, H is the mean curvature vector of the 
sphere in the spacetime, and dA is the volume element on the sphere. |8, 12 It has already been 



established that \^t,r \ = 47rr^, but it is also easily computed in this metric since 

dA = ^/\da^\ dO d(p = Vr^ sin^ 9 dO dip = sin 6 dO dip, (61) 

which yields that 

/• /'27r /'2n 

\^t,r\= / dA= / r^smededip= / 2r^ dip = ATrr^ . (62) 
Jj:t,r Jo Jo Jo 

To compute the mean curvature vector, H, we have that 

H = 7^'=n(5,-,9fc) (63) 
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where j, k & {9,ip}, j is the metric on E^ ,. (note that we have reused the variable 7 here to refer 
to the metric on the sphere and not to the metric on the t — constant hypersurface as we did in 
the section describing the framework of numerical relativity), and 11 is the second fundamental 
form tensor, which sends a pair of vectors tangent to the sphere to a vector normal to the sphere. 
It is defined as follows for all X,Y E TS^ 

JL{X,Y)^VxY- ^VxY (64) 

where V is the covariant derivative operator on N and is the induced covariant derivative 
operator on j.. Since 7 is diagonal, we only need to concern ourselves with the diagonal 
components of this tensor in order to compute H. To perform these computations, first note 
that the Christoffel symbols on the sphere in these coordinates have the following property (note 
that a superscript of 2 will always denote that that quantity corresponds to the sphere, also for 
the following computations, roman letters will denote subscripts in {9, </?}, while Greek letters 
will denote subscripts ranging over all four coordinates). 

^jk ~ lj9 ^ {.9jr],k~\~ Qkri,]^ 9jk,ri) = ~;j9 {,9jm,k~\~ 9km,j ~ 9 jk,m) = {,9jra,k~\~ 9km,j ~ 9jk,m) = Tj-fc 

(65) 

since both the radial and time directions are normal to the sphere. So we need not distinguish 
between the Christoffel symbols for the metric on the sphere and those for the entire manifold. 
We have then that 

n(9,, de) = Vede - ' V,9, = T,, ^ 9, - ^r,, ^ 9, = T,, * dt + T,, ' dr (66) 

and similarly 

u{d^, d^) = v^d^ - 'v^d^ = r^/' d, - dk = r^^* dt + r^/ a, (67) 

Then we need to compute the above Christoffel symbols. We have that 

^ee * = + 9er,,e - 9ee,n) = \9'\-9ee,t) = i^{-dt{r^)) = (68a) 

^ee = ^9'''^ {ger,,e + ger,fi - gee,r,) = 2g'''^{-ged,r) 

^-^^^^ = -r\\- — j (68b) 
r^^* = \9'''{9^n.^ + 9^n,^ - 9^^,n) = \9'\-9^^,t) = ^{-dt{r^ sin^ 6)) = (68c) 

^ - — ^ (-a,(r2 sin^ e)) = -r[\- — \ sin' Q. (68d) 



2 

Then we have that 



2M^ 






1 sin' 








2M 


r [\ - 






r 




2M 


r ( 1 - 






r 



n(a,, de) = r,, * dt + r,, ' Or = -r (^i - ^ J a, (69) 

n(5^, d^) = dt + r^/ dr = -r(l- sin' 9 9, (70) 
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This makes (63) become 



H = 7"''n(a9,c)„)+9wn(c)yS^, 



r 



1 



^2 sin^ 6 



r I 1 - — ] sin^ 9 dr 



(71) 



Finally, we can compute the Hawking mass as follows. 

1 



IGtt 



1 



167r ./s 



t,r 



^ H, H dA 



47rr2 / ^ 1 
IGtt 



IGtt 



r 



r 



r 



(9, 1 dA 



IGtt 



1 



JO 

27r PIT 



2M' 

1 ] g{dr, dr)r'^ sin 6' (i6' rfyj 



_2M 
2 V 47r V r 



r 
2M 



r 

27r /»7r 



sin 6* (i6' dy? 



JO 



sin 6* ^6* di^ 



t,r 



r f If 2M . 
2V ^4^1 ~ ~ I ^ 



r /2M 
2 V~ 
M 



(72) 

□ 



To prove Lemma 3.2, we will utilize the following additional lemma. 



Lemma A.l. For all rj G {t,r,6,(p}, let u-^j be as defined by equation (22). Then the following 
are true. 

1. For all 7] e {t, r, 9, ip}, Ric(z/^, z/^) = g{ur„ i^ri)g'''^ Ric^^. 

2. R = '^ g{Ur^, Ur,) Ric(z/^, Urj) 



Proof. Since {z/^} is a frame field, we have that 

Ric(z/t, ut) = e-2^ Ric(at, dt) = RiCfj = g{i^t, Vt)g^^ Ric« 
By similar arguments, we have that for rj G {t, r, 9, 

Ric(i/^, z/^) = 5((z/^, z/,,)^''''' Ric^^ 



(73) 



(74) 



18 



Alan R. Parry 



A Survey of Spherically Symmetric Spacetimes 



which proves [TJ To prove [2| we use the first result and find that since g is diagonal and 

g{uj^, z/^) = ±1, we have that 

R = Y.g'^^ Ric,^ = g^^ Ric,, = = ^ ^(^^, Ric(zy„ z/,). (75) 

Fact |2] also follows from a well understood more general property of frame fields and traces of 
tensors. □ 



With this, we prove Lemma 3.2[ 



Proof of Lemma 3_^. In this proof, all indices range over the coordinates {t, r, 6, By Lemma 



A.l and the well known formula for the components of the Riemann curvature tensor, R, in 



terms of the Christoffel symbols 

I ) m -p m -p m ip m p sp m 

^jkt — j£ ,k M ,j ~^ '- j£ ks k£ '- js ■ \'^) 

Since Ricjfc = R^-^^^, we have that 

Ric(z/j, ut) = g{i^u ^t)g^^ Ric« 

= e-'^ (r,,^, - r,,'^, + r,,T,/ - r,/r,/) (77a) 

_y I 2M 

Ric(z/f, z/r) = e Y 1 Ricfr- 



^' 1 ^ ij^tr'.k ^ kr ,t ~^ ^ tr'^^ ks ^ kr^ ts^) i^"^^) 



2M' 

Ric(l/r, I'r) = giyr) ^r)g^^ RiCrr = ( 1 ) RiCr 



(2M\ k k k k 

1 ~ j (^rr ,k " Tfcr ,r + ^ rr^ ks ~ ^ kr^ rs ) C^'^*^) 



Ric(z/e, i/e) = ^((z/e, z/e)fi'^^ Ricee = ^ Rice 



~ J.2 {^ee',k ^ke fi + Te/r^/ r^j/Fg/) (77d) 

Ric(z^<p, z/^) = 5f(z/<p, z^^)5''^^ Ric,.,. = ^ ^ ^ ^ Ric<p<p 

sm fc' 

"'■ (v k ^-lF ^f'^ — F *F ^\ (77p) 

ifiifi ,k kip ~ ipip ks kip ips ) • \' '^J 



"2 sin^ 9 



Then to compute these, we need the Christoffel symbols. For convenience, we will write the 
Christoffel symbols in four matrices, F*, F^, F^, and F''', which are the 4x4 matrices of Christoffel 
symbols corresponding to a given fixed upper index. To do so, recall, that the formula for the 
Christoffel symbols in terms of the metric components is the following 

^jk = -j^9'''"^{9jm,k + gkm,j - gjk,m) (78) 
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Using ( 78 ) and the fact that g is diagonal, we have that 

1 



r * 



jk 



9 {,9jm,k 9km, j 9jk,'i 



1 



9 {9jt,k + 9kt,j — 9jk,t 



1 



{9jt,k + 9kt,j — 9jk,t) 



2 5' {9jm,k + 9km,j 9jk,m) — ^9 ''{9jr,k + 9kr,j - 9jk,r) 



2M 



{9jr,k + 9kr,j 9jk,r) 



■ jk 



jk 



^9 "^{9jm,k + 9km,j — 9jk,m) — -9 {9je,k + 9ke,j — 9jk,e) 



^{9je,k + 9ke,j - 9jk,i 



(^9'^ {9jm,k + 9km,j 9jk,rn) — ^d'^ 

"^{9j'P,k + 9kip,j 9jk,ip) 



2r2 sin^ 6 



i9jtp,k ~l~ 9k(fi,j 9jk,ip) 



2r2 sin 6 



{9j'p,k + 9k^,j)- 



(79a) 



(79b) 



(79c) 



(79d) 



The last line is due to the fact that none of the metric components depend on ip. To help 
compute these quantities, it would be useful to note the following. 



9u,t = 9,(-e^^) = -2y,e^^ 
2M' 



9ttfi = <9e(-e 



2V\ 



9tt,r 

9rr,t = dtil 



9tt,i 





: 



(80a) 



9rr,e — de [ 1 — 



2M 



2M 



1 - 



dr l- 



2Mt 



2M 



2M 



-1 



1 - 



2M 



-2 



2M., 2M 



9ee,t = dt{r^)=0 

9ee,r = dr{r^) = 2r 

9^^,t = dt{r''^in^e) = Q 

9<fi<p,r = dr{r'^ sin^ 9) = 2r sin^ 9 



9ee,e 
9ee,ip 

9wfi 
9ipip,ip 



deir^) = 
d^{r^) = 

dg (r^ sin^ 9) = 2r^ sin 9 cos 9 
d^r^ sin^ ^ = 0. 



(80b) 



(80c) 



(80d) 



All the other metric components are 0. Using (79) and (80), it is a series of straightforward 
algebra computations to obtain the following. 





Vr 





o\ 
































^0 








0/ 



(81a) 
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1 - 



2M 



r 



1 - 





2M 



-1 



Mt / 2M 
r \ r 
2MYWMr M\ 
\ r r'^ ) 










r ( 1 - 



2M 







-r sin^ 6 ( 1 



2M\ 



/o 

















1 







1 


r 













r 








lo 





- 


sin cos 


/o 








1 











r 

cos 











sin 6' 




1 


COS 6 


I" 




r 


sin 9 



Using (81), we can now continue the computations started in (77). 

2M 



,2y 



+ Ke^^ 



r 



1 

2M 
r 



2-2V 



1 

2M 

r 

K- 
2M 



- dtVt - dt 



2M 



^ 2MY^ (Mr 
r J \ r 



M 



+ 



-,-2V 



^2V 



2M 



Mtt 



2M 



+ 2Vy 
2Mf 



1 - 



2M 



2M 

r 

2M 



r 



+ V'e'^ I 1 



- 2V:'e 



2J2V 



2M 
r 

2M 



1 - 



2V I Mr M 



+ 



+ Ke^^ 

VtMt 
r 

2Vre^ 
r 



2M 2Mr 



1 - 



1 - 



2M 

r 

2M 



r 



1 - 



2M 



-2 



(81b) 



^Ic) 



(81d) 
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Ric{ut, l^t) = (Vrr + V, 



+ 



2Vr 

r 

VtMt - Mu 



1 - 



2M 



+ 



Vr fM 



-Mr 



re 



2V 



1 - 



2M 



3Ml 

j.2q2V 



1 - 



2M 



(82) 



^ \l^-—{dtiVr)+dr\^[l-—j ) -dt{Vr) 

dtlil 1 { — - — ]\ -2dti-]+Vr{Vt + 



_y 2M \ ^ ^ f 1, 2M\\ ^^f I, 2M 



= e-^ y 1 - — ja.a, log (^1 - — J ) - ^t^^ ( -o log ( 1 
2Mt ( 2M\~^ 



2Mf ( 2M^~^/^ 



^<^t.Vr)^~:^[\-~-^] (83) 



1 - 



2M 



dt 



^2V 



1 - 



2M 



-drVr-dr 1- 

2M 



2M\~WMr M 



+ 



+ 1 



1 

2M 

r 

r 



— - 



1 - 



-2a 



2M 



-1 



+ 



/Mr M 
\ r 



2y 



1 - 



2M 



+ Vr + 



-2 / 2 



\ r r'^ ) 

'M,. M^ ^ 



^ 2M\~^ (Mr M\ 
r J \ r J 
2M2 ' 



y^2^2 



1 



2M 
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2M\~^ / 2M\' 

r J f'^Q^v y J. J 



2 Vr\ fM 



VtMt - M, 



} 



I \ \ Til \ sm y / \ T I \ r 



[Mr M\ ( 2M\ cos^e 
-r — +2 1- 



\ r r'^ J \ r J sin 6 
-i -fl-— Vrf— -^V — -rfl-— 1 fv; + ^^ 

'M Mr\ ( 2M\ cos'^e 
+ r^-— +2 1- ^ 



r y \ r J sin 6* 

1 I sin^^ 2M M 1 f 2M 

^{ - + + + - ri; I 1 - 

""^ ' sm y r r sm 

2M\ / 2M\ cos'^e 
-2 1 +2 1- 



sin^ ^ I \ V ^ / 



Or -rsin'e 1 +ae(-sin^cos^) 



r / sin ^ 

1 /M \ Vr f 2M\ , , 

Ric(^,, ve) = y\-^^r)-^[y-—) (85) 

1 \ ^ ( 2/1 A 2M^^ 



- r sm^ e[l \Vr + - \ - r sm^ 0[— -] - cos^ 

r J \ r J \ r 

+ 2 sin^ ^ ( 1 - — 1 + 2 cos^ Q 



23 



Alan R. Pcirry 



A Survey of Spherically Symmetric Spacetimes 



sin^ eil r sin^ Oi^ - cos^ 9 + sin^ 9 



sin^ ^1 \ r J \ r 
-rsin^^ ( 1 1 ( K + - 1 +rsm^e[ — -] - cos^ 9 



+ 2 sin^ e{l- + 2 cos^ 



2M\ / 2M\ ( 2M\ 
-21 +2 1 



+ M, -rKU -21 +2|i I 

Ric(i.^, i.^) = ^ (^y + - 7 (l - ^) = ^e). (86) 

These are the only nonzero components of the Ricci curvature tensor since 

1 „. 1 



1 



j.qV 



te ,k k9 ,t ~^ te ks ^ kO ^ts ) 



^ ( p r p 6 p p rp A; p Sp fe p V'p k\ 

j.qV \ rd ,t ee ,t ipe ,t '- ke '- tr ^ ko ^te ~ ^ ko ^ ) 

/ p rp 9 p 9p r p V'p 



(87) 



Ric(i/t, i/^) = ■ ^ y Rict^ = ■ ^ T/ R-tfcy 



rsm (/^e" 

/p fc p fe _|_ p sp p sp 

ti^ ,fe kif> ,t ^ tip ks ktp ts ) 
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r sin 6'e^ V ^ w ^ ipip ^ td ^ np ^ tip ^ Oip ^ tip J 
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/ 


2M 








r 


1 


2M 









r 



(p /c P ^ _1_ P ^P ^ P ^P 

re ,k k9 ,r ~^ re ks ^ fe6» ^ rs j 

/p 6 P ''^ _L P ^P P ^P ^\ 

V re ,9 ~ k9 ,r r9 ke ~ k9 rs ) 



pr p'^ P'^_lP ^P _l P ^P ^ 
r6> ,r ~ ^ ,r ~ ^ (^0 ,r ^ r6> ^ r6l + ^ rS ^ 6>0 



I p 0p p sp r p sp 9 p sp 

^ r6» ^ re '-rs ^ ee ^ rs ^ ipe ^ rs ) 
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1 /, 2M /ms^_ , , 



11^ 2M f cos 9 
r " 




cos 6 



r \r sin 9 r sin 9 



(89) 



r sin^ 
1 

r sin 6* 
1 

r sin 6* 
1 

r sin6' 



2M 
1 Ric 



r</3 



r sin ^ 



2M, 



R 



'rkip 



1 - 



2M 



r 



/"p /i: "p fc _i_p ^p ^ p -^p ^^ 

V ri^ ,fc fei/p ,r "I" rip fcs fcip rs ) 



i-^(-r ^ -r 



r 



^pip ' vi^ Tip ' vi^ Oip 



I p (^p i/:! p sp r p sp p sp 

Tip ^ip rs Oip rs ip)ip 



Lpip rs 



r sin9 




1 - 



2M 



( —Y ^' — F ^ — F "^Y — Y ^Y "^1 



(90) 



Ric(z/e,i/^) 



1 



Ricfl 



^ T> k 

7^9 



/ p "P ^ _L "P "^"P ^ P "^P ^\ 

K'- Oip ,k ^ kip ,e ~^ ^ e<p ^ ks ^ kip ^ es ) 

/p ifi p k _|_ p Vp p sp fc\ 

l.-"- 693 ,ip ktp ,e ~^ 8ip ktp ktp es ) 



sin 9 
1 

sin ^ 
1 

sin 9 

^ f -r »• _p ^_p "^-lf "^f ^ _l f "^f 
r2 sin 6* ^ '"'^ ''^ J- (^^^ ,e + J- '"v' + ^ 6 

I p i/'p p ''p p ^p p v'p 

"T 6l(/3 kip Or kip ee kip Op ) 

( p rp 9 p 0p r p i^p (ysN 



sin 6' 




(91) 



The other components are by the symmetry of the Ricci curvature tensor. To get the last 



statement, we note that by the last line above and Lemma A.l 2, we have that 

R = g{i^ri, ^T]) Ric(z/^, z/^) = — Ric(z/f, Ut) + Ric(z/r, Ur) + 2 Ric(z/e, t^e)- (92) 



From here, it is a matter of algebra and the statements above to get the desired equation for 
R. □ 



Next, we prove Proposition 3.7 
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Proof of Proposition 3.1. We already have the equations resulting from two of the components 



of the Einstein equation, namely (30) and (36), which we reprint here for convenience. 



Vr 



1 - 



2M 



-1 



M 



+ AirrP 



(93) 
(94) 



In order to show that solving these two equations is sufficient to solve the entire Einstein equa- 
tion, we must show that if these equations hold, so do the other components of the Einstein 
equation. To do so, we ffist need to write down the other components of the Einstein equation. 
Recall that the above equations come from the ut, vt and z/^, Vr components of the Einstein equa- 
tion. Thus we only have the vt, and z/^, vq components left to compute (the z/<^, v^p component 
is identical to i/g, vq component). Then, by Corollary 3.3 and (42), we have the following. 



2Mt 



2M 



-1/2 



inp 



2M 



(95) 



and 



8ttT{u0, Ug) 



Vrr + K + 



+ 



r 
Mu 



1 - 



2M 



re 



2V 



r 

2M 

r 



+ 



M 
r 

3Ml 

^2 g2 



Mr 



1 Vr 



2M\' 
r J 



(96) 



Next, since T is a stress-energy tensor, we know that it is divergence free, that is, Vg ■ T = 0, 
where Vg-T denotes the divergence of the tensor T with respect to the metric g. We will compute 
next what this equation means in terms of the functions /i, P, p, and Q. In what follows, 
the summing indices run through the set {t,r,6,(f} and we will use the Einstein summation 
convention wherever it applies, while specifically denoting any other summations of a different 
form. First, we define 

for all k G {t, r, 6, (f}. Then we can write the divergence of T as follows. 



(97) 



iVu,T){ut, d,) dx' + {V.J){vr, dj) dx' + (V.,T)(z/e, dj) dx' + (V.^T)(z/^, d,) dx' 

(9^ 



We will simplify each of these four terms individually. By equations (81), (22), and (42), we 
have that 



(V,,T)(z/„ dk) dx' = e-2^(ViT)(ai, d^) dx 



-2V 



dtindt, d,)) - T,rndm, a,) - r,,'"T(a„ 



dx 
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= e 



-2V 



dt{T{dt, dt)) dt + dt{T{dt, dr)) dr - 2T,,'T{dt, dt) dt - 2T,/T{dt, dr) dt 
(r,,* + V,;)T{dt, dr) dr - r,^'T{dt, dt) dr - r,/T{dr, dr) dr 



— e 



-2V 



2pVre^^ 



. 2M , 
1 dt 



2M 
r 

M, 



dr - 2Vte'^^n dt 



14 + ^(1-^ 



eV 1 



2M 



-1/2 



dr 



VrQ^^ndr-Vre'Tdr 



2V , 



— e 



-2V 



2M 



dt 



+ Vte^'p 1 



2M\ 
r J 



-1/2 



+ e^A ( 1 



+ 



r \ r 



eV 1 



2M 



-1/2 



Vre^^p - Vre^^P 1 dr 



2M 



2V 



-3/2 



2A/f \ / / 2A/f^ ~^/'^ 

(V.,r)Ka,)dx^= (/i,-2pKe^-^/l-— dt+ p,e-^h - — 



(99) 



Next, we have 
{V.^T){iyr,dk)dx' 



1 - 



= 1 



2M' 
r 

2M^ 
r 

2M^ 



{VrT){dr,dk)dx'' 
dr{T{dr, dk)) - T^;^T{dm, dk) - r^k^'Tidr, dm) 



dx^ 



-2V * 

rr 



dr{T{dr, dt)) dt + dr{T{dr, dr)) dr - (F,/ + T,/) T{dt, dr) dt 

'T{dt, dr) dr - 2V,;T{dr, dr) dr - VjT{dt, dt) dt - r,/T(9„ dr) dt^ 



1 - 



2M~ 



a, pe^ 1 1 \ dt^dr\P\l- 



-^ 



dr 



-peMl- 



2M\"^/^ 



1 - 



2M 



2pMt A _ 2M\ 



re^ 



-5/2 



1 - 



^ dr - 2P (^1 - — - 

2M\"^ , PMt /. 2M 



2M 



dt 



Mr M . , 



1 - 



dt 
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= 1 



2M 



2M 



-3/2 



Mr M 



+ pVre' 1 



2M^ 



-1/2 



+ p.e'' 1 - 



PMt 2M 



Pr\l 

2pMt 



2M 



2M 



-1/2 



1 



r \ r / r 

2M\~WMr M 



Vr 



dt 



2M 



re" 



+ 2P ( 1 

2M\"^/^ / 2M 

1 - — j - 2P I 1 



V /-, 2M Mi 



2pMt 2M 



re 



V 



2M 
r 

-3/2 



(p + P)] dt 



dr 



Thirdly, 



{V,,T){ue,dk)dx^ = -jVeT){de,dk)dx' 



dx^ 



1 



r,/T(9„ 9t) dt - T,;T{dr, dr) dr - r jT(a,, de) dr 



rpe 



2M 



dt + rP dr — rQ dr 



{V.,T){ue,dk)dx' = ^ 



Lastly, we have that 
{V,T){iy^,d,)dx'^ 



^ I 2M , P-Q , 
1 dt+ dr 





• 2 /) V 

sm 




1 




sin^^ 




1 




sin^^ 




_ r f 








1 



iW^T){d^,dk)dx'' 



(100) 



(101) 



dx'' 



- T^;T{dr, dt) dt - V^;T{dr, dr) dr - VjT{de, dg) de 



sin^ 9 



^ifj dip) 



2M 



rpe sin 9\j 1 dt + r sin 9P dr 



+ r'^Q sin 9 cos 9 d9 - rQ sin^ 9 dr - r'^Q cos 9 sin 9 d9^ 
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dt H ^ dr 



{V,T){u^,dk)dx' = {V,,T){u0,dk)dx 
Then equation ([98|) becomes 



(102) 



V, ■ T = -{V,J){ut, dj) dx^ + {V,T){ur, d,) dx^ + {V,,T){ue, dj) dx^ + {V ,T){v^, d,) dx^ 

-1/2 



2M 



dt 



2M Mt 



+ Pr 



2pMt 



re' 



-Pt + ^pe^xn 



2M 



1 - 



-3/2^ 



2M 



r 



2M 



{p + P)\ dt 



Vr{p + P) dr 



dr+'-^Jl-'-^dt+'-^^dr 



-V 



2M 



2M 



Vr + - ] + Pre 



2M Mt 



-1/2 



Vrip + P)+Pr 



2pMt 



1 - 



r 
2M 



r 



re 



V 



2M 



{p + P)j dt 

2{P-Q) 
r 

(103) 



dr 



Then V„ ■ T = yields the following two equations 



pt = 2pe \/ 1 - 



2M 



Vr 



pre 



V 



1 - 



2M Mi 



1 - 



2M 



pt 



2M 



Vr{p + P)+Pr + 



2(p-g) 



r 

2pMf 
r 



2M 



iP + P) 

1 



(104a) 
(104b) 



We are now ready to show that solving (30) and (36) on each t = constant slice also solves 
the Einstein equation. In other words, as stated earlier, we want to show that solving these 
two equations for each constant t automatically satisfies the remaining two unique nonzero 
components of the Einstein equation, (95) and (96). We will show first that (30), (36), and 
(104) implies that (95) holds as well. To do this, we first need to show that (30) and (95) are 
compatible, that is, that, given these equations, Mrt = Mtr- Differentiating (30) with respect to 
t yields. 



dtMr = dt [A-nr'^p] 
Mrt = A-nr^Pt 



(105) 
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while differentiating (95) with respect to r yields, 



drMt = dr 47rr^e^pWl - 



2M 



Mtr = 47r 



2re^p\l 1 - — + r^Vre^pd 1 - — 



+ r e 1 r e p i — 

r \ r J \ r 



M 



To make the simplification process easier, to the right hand side above, we will add in and sub- 
tract out a copy of the second term in the sum above, namely, the term Aur'^Vre^ p^Jl — 2M/r. 



tr 



An 



2re>Wl 



r 



2M , ^ / 2M . V , 
2rV,e^pWl rV.e^pWl 



r 



2M 



+ r'^e^pr 



2pe 



Mt 
/it H 



V . 





2M 








r 




2M 








r 



r^^'p I 1 



2M 



(Mr 
\ r 



M 



rp I I ~|~ 

r 



V 



Airr 



2M 



r 



e^pll 



-1 



2M 



r 



1 

1/2 



2M 



Mr M 



2M\ , 
1 1 (p + P) 



e^pll 



2M 



M, 



1 - 



r 

-1/2 

2M 



M 



+ AnrP 



e^pll 



(p + P) - 47rrpe^ ( 1 - 



2M 



-1/2 



r 

2M 
r 



^47rrp 



M 



-1/2 



(/i + P) 



47rr pt 



(106) 



where in the next to last line we made substitutions using equations (30), (36), and (104a) and 
the last line used (95). Thus if ■ T = 0, (30) and (95) are compatible, that is, Mtr = Mrt 



everywhere the equations are defined, that is, wherever r 7^ 0. However, by using L'Hopital's 
rule on the equations above coupled with the fact that at r = 0, M = Mr = Mrr = for all 
t, the equation Mtr = still holds at the central value r = 0. This implies that there exists 



a function M{t,r) which satisfies both (30) and (36) everywhere, which of course is the metric 
function we seek. 



If, for all values of t, we integrate (30) with the initial condition M = at r = 0, we 



will obtain some function M*{t,r) which satisfies (30) everywhere. Then we have that M, 



Mr everywhere, where M is the function that satisfies both (30) and (36). This implies that 



M*{t, r) = M{t, r) + f{t) for some smooth function f{t). However, we also have that M*{t, 0) 
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M(t,0) = for all t, which implies that f{t) = 0. Hence M*{t,r) = M{t,r). Thus the function 



we obtain by integrating (30) with compatible initial conditions necessarily also satisfies (95). 



Finally, we use equations (30), (36), (104), and (95) (since the first three imply it) to show 



that (96) is automatically satisfied. In order to do this, we will have to compute Vrr and Mu- 



They are 



dtMt = dt I ATcr^e^'pdl 



2M 



V^te^PA/1 



2M 



+ eVtA 1 



2M e^pMt 



1 - 



2M 



2M pMt f 2M 
pt\ 1 ( 1 



-1/2- 



'1/2' 



(107) 



and 



2 1 



1 - 



2M 



-1 



M 



+ AnrP 



2M 



2 



/Mr M 

\ r 



M 



+ A-KrP I + ( 1 



2M 



-1 



Mr 2M 



2MYWMr M\ , 

2K ( 1 — - ^ + I 1 

r J \ r r-^ ) 

2M\'^ (Mr M\ Vr . 
2K(1 -^+|1 



2MY^ (Mr M M 



+ 47rP + 47rrP. 



^3 ^3 



+ 47rP + 47rrP,. 



2M 



Mr M 



Vr 



+ 



1 



2M 



Mr 

r / 



M\ Vr 



^ + 47r 1 



2M 



+ SvrP + A-nrPr 



(2P + rPr 



Making these substitutions into (96) yields 

Vr 



Vrr + V' + 



+ 



VMt - Mt, 



re 



2V 



1 - 



2M 



2M 



M 
r 

3Ml 



Mr 
1 - 









r 


2M\ 


-2 



4.(2P+rP.) + K=(l-?^) + 'ifM„-^) 

r / r \ r I 



A 1- 



2M\ pMt 2M 



r 



r 



r 



Atx{2P + rPr) + Vr\^+ AnrP 



M 



Ptil 



2MY^''^ pMt 



2M 



-3/2 



M\ 
r J 

-3/2 



3Ml 



UnpMt 



2M 



2M 



'3/2 
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4vr(2P + rP.)+4vrrKP+^-'^ 

r 



47r(2P + rPr) + 47rrV;(P + fi) 
47r(2P + rPr) + 47rrV;(P + fi) 

SttQ 



Anrpt 
871 pMt 



1 - 



2M 



2M 



^1/2 



2M' 



-3/2 



2M 



2{P-Q)\ SnpMt 



'1/2 



-3/2 



SnpMt 



1 - 



2M 



-3/2 



2M 



-3/2 



(109) 



where the last two lines follow from using ( 104b ) . Thus so long as equation ( 104 ) holds, equations 
(30), (36), and (95) imply equation (96). Since equations (104), (30), and (36) imply equation 
(95) and equation (104) follows from ■ T = 0, we have that Vg ■ T = implies that solving 
equations (30) and (36) at each time t is sufficient to solve the entire Einstein equation, which 
was the desired result. □ 



The proof of Lemma 4.1 follows. 



Proof of Lemma \4.1\ Recall equation (47). Then to compute the components of T, we ffist have 
that 



df = ftdt + fr dr = pe \ 1 



2M 



df = ftdt + fr dr = pe^\l 1 
Then, since g is diagonal, this implies that 



2M 



dt + fr dr 



dt + fr dr. 



(110) 
(111) 



M/r = 9{df, df) = 9 { ^/ 1 - ^ rft + dr, pe^ \l 1 



2M 



dt + fr dr 



2V 



2M 



ppg{dt, dt) + frfr g{dr, dr) 



2M 



ilfrf-lpf). 



(112) 



With these facts, we now compute the quantities in question and do so in the same order as 
they were presented. Thus we have the following. 



T{ut,^t) = /io I Y2^f{'^t)df{iyt) 



1 



2M\ \frf-\pf 



■^2 



+ i/r 



g{i^t, ^t) 



2 

■^2 



2M 



PP+I/I + 1- 



2M\ \fr\'-\p[ 



pA\ft+{i- 



2M\ \ fr\' + \pf 
f2 



(113) 
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Next, because Ujj is an everywhere orthonormal basis 
1 



/^ol ^1 1- ^ 



{Pfr + frP) 



(i - ^ ) Mm 



2 / 2M\ 
/iO I 1 Jr/r 



2M\ iM'-br 



■^2 



+ i/r 



^2 



r 



■^2 



2M\ \frt + \p[ 



(114) 



Finally, 



^2 



-/io |/|'+ 1 



2M^ I M^-IpI^ ' 

^2 



^2 



r 

2 I |2 



so that T(z/<^, z/^) = T(z/6), 



Finally, we prove Lemma 4.4 



Proof of Lemma \4-4\ By equation (42) and Lemma 4.1, we have that 

2M\ \frf + \pf 



(115) 



/i = /io i/r+ 1 



■^2 



^ = /^o -i/r 



Q = -/io 1 



2M\ iM' + br 



^2 



2M^ \ frf-\pf ' 
■^2 



(116) 

□ 



(117a) 
(117b) 
(117c) 

(117d) 
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We know from a previous proof that ■ T = is equivalent to equation (104). Moreover, by 
(48) and Lemma 4.3, we have that (46) is equivalent to (54) and (58). Thus it suffices to show 
that given (117), equations (54) and (58) imply equation (104). 



We will start with (104a). On one hand, the left hand side is equal to differentiating (117a) 
with respect to t. 



l^t = dt 



^0 



/io |/|'+ 1- 



2M\ \ frf + \pf ' 

^2 



ftf + fft 



2Mt i\fX + \pf 



^2 



1 - 



2M\ frtfr + frfrt + PtP + PPt 



e^/l-^(p/ + /p)- 



2Mt (\f,f + \pf' 



^2 



-I fV h 2M Pi 



2Mi (\f,f + \pf 



1 - 



2M 



^2 



2M\ frtfr + frfrt + PtP + PPt 

2M 



f2 



2M 



/,aJpe^./l 



2M 



+ frdr(pe''Jl- 



2M 



:il8a) 



On the other hand, we can substitute (117) into (104a). 



fit = 2pe^ W 1 - — ( K + - 1 + Pre' \l 1 



2M Mt 



2M 



ip + P) 



2e'\ 1 



2M / 2M 



■^2 



r V T2 

2M 



1 - 



MfrP) ][Vr + 



r 



^e-(l- 



-f- r^2^r 



1 

2M\ 



r 

3/2 



Re{frp) 



, 2M 2Mtfio ( \fr\' + \pf 



r 



r 



2M\ 



ifrP + pfr)[Vr + 



r 



ifrP + Pfr) 



/l_2M 2Mtflo (\frf + \pf 



r 



r 



■^2 



Pf2 



+ P 



2M 



frVre' 



2M 2fre 



V 



+ 



1 



2M 



■^2 



2M 



fr^r^ 



ifrP + Pfr)dr f 1 



V 



V 



2M 2fre 



+ 



r 



1 



2M 



2MtfIo ( \fr\'+\pf 



■^2 



2M 



r 

2M 

r 



^2 

2M 



i,frrP ~l~ Pfrr ~l~ frPr ~l~ Prfr 
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_/io / 2M 



frVre\ l 



2M 



2M 2fre 



V 



+ 



1 



2M 



2M 



-|- frr(i 



V 



V 



2M 2fre 



+ 



1 



2M 



r 



r 



■^2 



2M 



-1/2 



9. 1 



2M 



2M 



2M 2fre 



2M 



+ e \1 ifrPr+Prfr) 



+ 



1 



2M 



fio A 2M 



2M 



+ frrG \ l 



V 



2M 2fre 



+ 



2M 



Ke^A/l-^(/,p + p/,) 



\ 'T' 

2M,/io (\frf + \pf 



2e\frP + Pf r)dr \1 



2M 



r 



2M 



ifrPr +Prfr) 



/it = Pj^ I 1 



2M 



9. l/.e^Wl 



2M \ 2fre 



V 



2M 



/xo A 2M 



2M \ 2fre 



V 



+ 



r 



2M 



Y'2 



frdr\pe^\ll- — ]+frdr\pe^\ l 



(118b) 



Setting ( |118a[ ) equal to ( |118b[ ), we obtain 

2M^ 



1^0 



P\fe''\ll- 



2M ^ 



+ p fe\l- 



r 



2M 



2M, (\frf + \p\' 



Y2 



_fio ( 2M 



/re \/l 



2M \ 2/^e 



V 



+ 



1 



2M 



/io A 2M 



5J/.e^A/l- 



2M \ 2/^e 



V 



+ 



2M 



■^2 



^9.|pe^/l-^ +^5Jpe^/l-^ 
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P I TVe^ ( 1 - 



2M 



-1/2 



2M 



-1/2 



+ Pt 



V 



2M \ 2fre^ 



2M 



P 



dr fr-e \ 1 



2M \ 2fre 



+ 



2M 



= p 



+ p 



2M 



-Pt + -T^f 1 - 



-1/2 



+ 



Vr 



r 



2M 



r 



dr frS \ l- 



2M 



2M 



r 



r 



:il8cl 



Note that, since V and M are real valued, the second term above is the complex conjugate of 



the first term. Then by (54) equation (118c) holds. 



Next, we consider (104b). We first differentiate (117b) with respect to t and then compare 



it to ( |104b[ ) after substituting ( |117[ ), ( [50| ), and ( |5l| . This yields, 

2/io 



Pt = dt 



1 - — 1 Reifrp) 



r 



f2 



4M 



r 

4Mt 
r 



Re(/^p) 



Re(/,p) + 1 



2M 
r 

2M 

r 



{frtP + frtP + frPt + frPt) 



pdr\pe \ I- 



2M 



+ pdr ^pe \/l 
/io I 4Mt 



2M 



^2 



Re(/,p) + 1 



+ frPt + frPt 
2M 



pdr \pe \ 1 



2M 



+ Prpe^\ll - ^ + \pfdr I e^^/1 



2M 



+ frPt + frPt 



■Y-2 



4Mt , / 2M 
-^Re{frp)+{l 



2M 



+ \p\'dr\e\l 



i2 o / V I - 2M 



/rPt + frPt 



(119a) 
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From (104b), we obtain 



2M 



Vrifl + P)+Pr 



2{P-Q)\ 2pM, 



1 - 



2M 



1 



2M 



2/ioK 

■^2 



1 - 



2M 



i\frf+\pf) 



+ fiodr -|/|'+ 1 



2M\ \f/ + \pf\ ^ 4/xol/. 



■^2 



+ 



rT2 



2M 



1 - 



2M 



2M 



2KIM \/i + 2Kbr\/i 



2M 



2M 



-1/2 



+ 1 



2M 



-1/2 



a 1- 



2M 



+ 



1 - 



2M 



2M 



y-r2 . 



2M 



-1/2 OoV. 



y r 



+ 2Ke^ W 1 - + 2 e^9, I W 1 - 



2M 



2 



2M 



+ 1- 



2M 



+ 2|M^^a,, I ^/l 

2M 



^-2 



1 - 



2M 



<9r /re^A/ 1 



2M 



+ e^Wl (frrfr + frfrr) }- 



rT2 



Re(/^p) 



2M 



+ eM-TV 1 



2M 



-1/2 



^2^ / _2M 



dr fre \ I 



+ fr 



2M 



2M 



-1/2 



2M 



5. I H'e^Wl 



2M 



4/ioMt 



Re(/^p) 



(119b) 
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Equating (119a) and (119b) yields 

2M" 



/io / 2M 

^2 ' 



5rl N'e^Wl- 



fr 



2M 



r 

-1/2 



2M" 



+ -T^/ 1 



a. /.e% 1 



2M 



+ e 



-t7 1- 



2M 



+ \p\^dr\e\l 



,2c / y 2M 



2M 



r 



2fr 



1 



2M 



-1/2 



^2^ / _2M 



4/ioM, 
rT2 



Re(/^p) 



= /. 



2M 



-1/2 



2U L_2M 



2M 



r 



+ fr 



-Pt + e 



-TV 1 



2M 



r 



r 



r 



:ii9c) 



As before, since M and are real valued, the second term above is the complex conjugate of 



the first term. Then by equation (54), equation (119c) holds. Then equations (54) and (58) 



imply equation (104) holds, which completes the proof. 



□ 
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